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Abstract 

We describe recent theoretical and experimental progress on mak- 
ing objects invisible. Ideas for devices that would have once seemed 
fanciful may now be at least approximately realized physically, using 
a new class of artificially structured materials, metamaterials. The 
equations that govern a variety of wave phenomena, including elec- 
trostatics, electromagnetism, acoustics and quantum mechanics, have 
transformation laws under changes of variables which allow one to 
design material parameters that steer waves around a hidden region, 
returning them to their original path on the far side. Not only are 
observers unaware of the contents of the hidden region, they are not 
even aware that something is being hidden; the object, which casts no 
shadow, is said to be cloaked. Proposals for, and even experimental 
implementations of, such cloaking devices have received the most at- 
tention, but other devices having striking effects on wave propagation, 
unseen in nature, are also possible. These designs are initially based 
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on the transformation laws of the relevant PDEs, but due to the sin- 
gular transformations needed for the desired effects, care needs to be 
taken in formulating and analyzing physically meaningful solutions. 
We recount the recent history of the subject and discuss some of the 
mathematical and physical issues involved. 



1 Introduction 

Invisibility has been a subject of human fascination for millennia, from the 
Greek legend of Perseus versus Medusa to the more recent The Invisible Man 
and Harry Potter. Over the years, there have been occasional scientific pre- 
scriptions for invisibility in various settings, e.g., [56j [11]. However, since 
2005 there has been a wave of serious theoretical proposals [U E3J UJE E3 E3] 
in the physics literature, and a widely reported experiment by Schurig et 
al. [99], for cloaking devices - structures that would not only make an ob- 
ject invisible but also undetectable to electromagnetic waves, thus making it 
cloaked. The particular route to cloaking that has received the most atten- 
tion is that of transformation optics |116j . the design of optical devices with 
customized effects on wave propagation, made possible by taking advantage 
of the transformation rules for the material properties of optics: the index 
of refraction n(x) for scalar optics, governed by the Helmholtz equation, and 
the electrical permittivity e(x) and magnetic permeability for vector 
optics, as described by Maxwell's equations. It is this approach to cloaking, 
and other novel effects on wave propagation, that we will examine here. 

As it happens, two papers appeared in the same 2006 issue of Science with 
transformation optics-based proposals for cloaking. Leonhardt [73] gave a 
description, based on conformal mapping, of inhomogeneous indices of re- 
fraction n in two dimensions that would cause light rays to go around a 
region and emerge on the other side as if they had passed through empty 
space (for which n(x) = 1). On the other hand, Pendry, Schurig and Smith 
P3] gave a prescription for values of e and \i yielding a cloaking device for 
electromagnetic waves, based on the fact that e and fi transform nicely un- 



der changes of variables, cf. (25). In fact, this construction used the same 
singular transformation ^ as had been used three years earlier [M| H5] to 
describe examples of nondetectability in the context of Calderon's Problem 
for conductivity, which transforms in the same way as e and [i. 
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We briefly outline here the basic ideas of transformation optics, in the context 
of electrostatics, leading to a theoretical blueprint of a conductivity that 
cloaks an object from observation using electrostatic measurements [HI l4"5] . 
Given that the invariance of the underlying equation is a crucial ingredient 
of transformation optics it is natural to set Calderon's problem on a compact 
Riemannian manifold with boundary, (M, g) with g the Riemannian metric 
and boundary DM where the observations are made. The Laplace-Beltrami 
operator associated to g is given in local coordinates by 



where (g^) is the matrix inverse of the metric tensor (g^) and \g\ = detg. 
Let us consider the Dirichlet problem associated to 0, 

A g u = on M, u\ dM = f. (2) 

We define the Dirichlet-to-Neumann (DN) map in this case by 



i,7=l v J/ 



(3) 

8M 



where v denotes the unit outer normal. Calderon's (inverse) problem, the 
question of whether one can recover g from A g , has been the subject of a 
tremendous amount of work over the last quarter century. In Sec. 2, we 
briefly summarize the history and current status of this problem. 

Given the invariant formulation of the DN map, it is straightforward to see 
that 

K V g = Ag (4) 

for any C°° diffeomorphism ip of M which is the identity on the boundary. As 
usual, ijj*g denotes the pull back of the metric g by the diffeomorphism ip. For 
domains in Euclidean space of dimension n > 3, the metric g corresponds to 
an anisotropic conductivity a, represented by the symmetric matrix-valued 
function 

= \g\ 1/2 g ij . (5) 
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The DN map sends the voltage potential at the boundary to the induced 
current flux. 

The invariance Q can be considered as a weak form of invisibility. How- 
ever, although the (generally distinct) media represented by g and ip*g are 
indistinguishable by boundary observations, nothing has yet been hidden. In 
cloaking, we are looking for a way to hide from boundary measurements both 
an object enclosed in some domain D and the fact that it is being hidden. 
Suppose now that an object we want to cloak is enclosed in the ball of radius 
one, -8(0, 1), and that we measure the DN map on the boundary of the the 
ball of radius two, .6(0,2). Motivated by degenerations of singular Rieman- 
nian manifolds (see Sec. 3) consider the following singular transformation 
stretching (or "blowing up") the origin to the ball -6(0, 1): 




Also note that the metric g = (Fi)*g , where (-F\)* = (F^ 1 )* and g is the 
Euclidean metric, is singular on the unit sphere § n_1 , the interface between 
the cloaked and uncloaked regions, which we call the cloaking surface. In 
fact, the conductivity a associated to this metric by (|5| has zero and/or 
infinite eigenvalues (depending on the dimension) as r \ 1. In IR 3 , a is given 
in spherical coordinates (r, 0, 9) \— > (r sin 9 cos 0, r sin 9 sin 0, r cos 9) by 
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/ 2(r - l) 2 sinfl \ 

a = 2sinfl , 1 < \x\ < 2. (7) 

\ 2(sin#)- 1 / 

Note that a is singular (degenerate) on the sphere of radius 1 in the sense 
that it is not bounded from below by any positive multiple of the identity 
matrix I. (See [62] for a similar calculation.) 




Figure 2: Analytic solutions for the currents 



The currents associated to this singular conductivity on 5(0,2) \ 5(0,1) 
are shown in Fig. 2. No currents originating at (95(0,2) have access to 
the region 5(0,1), so that (heuristically) if the conductivity is changed in 
5(0,1), the measurements on the boundary 95(0,2) do not change. Any 
object in 5(0, 1) is both unaffected and undetectable by currents from the 
outside. Moreover, all voltage-to-current measurements made on 95(0,2) 
give the same results as the measurements on the surface of a ball filled with 
homogeneous, isotropic material. The object is said to be cloaked, and the 
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structure on 5(0,2) \ 5(0, 1) producing this effect is said to be a cloaking 
device. 

However, this intuition needs to be supported by rigorous analysis of the 
solutions on the entire region -8(0, 2). If we consider a singular metric g de- 
fined by (.Fi)*^) on 5(0, 2) \5(0, 1), an arbitrary positive-definite symmet- 
ric metric on 5(0, 1), and 5 X (5(0, 2)) smooth solutions of the conductivity 
equation, it was shown in |HJ HH] that, for n > 3, the following theorem 
holds. 

Theorem 1.1 = A go . 

In other words the boundary observations for the singular metric g are the 
same as the boundary observations for the Euclidean metric; thus, any object 
in B(0, 1) is invisible to electrostatic measurements. We remark here that 
the measurements of the DN map or "near field" are equivalent to scattering 
or "far field" information [12]. Also, see [62] for the planar case, n = 2. 



In the proof of Thm. L_l one has to pay special attention to what is meant by 
a solution of the Laplace-Beltrami equation (J5]) with singular coefficients. In 
[HI US] , we considered functions that are bounded and in the Sobolev space 
5^5(0, 2)), and are solutions in the sense of distributions. Later, we will 
also consider more general solutions. 



The proof of Theorem 1.1 has two ingredients, which are also the main ideas 



behind transformation optics: 

• The invariance of the equation under transformations, i.e., identity Q. 

• A (quite standard) removable singularities theorem: points are remov- 
able singularities of bounded harmonic functions. 

The second point implies that bounded solutions of the Laplace-Beltrami 
equation with the singular metric indicated above on the annulus 5(0, 2) \ 
5(0, 1) are equivalent to bounded harmonic functions on the whole ball 
5(0,2). This shows that any H 1 solution u to the equation ([!]) is con- 
stant on the ball of radius 1 with the constant the value of the corresponding 
harmonic function v (0) with v — u o F\. 

The 2003 papers |331SH] were intended to give counterexamples to uniqueness 
in Calderon's problem when the anisotropic conductivity is allowed to be 
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only positive semi-definite. In the summer of 2006, Bob Kohn called our 
attention to the paper [93] where the same transformation F\ was used to 
propose cloaking for Maxwell's equations, justified by the analogue of (El). 
In fact the electrical permittivity and magnetic permeability in the blueprint 
for a cloaking device given in [03] are 



(&) = m = {\g\ x/2 g ij ) (8) 

with g = (Fi)*g . The proposal of [73] (appearing in the same issue of Sci- 
ence]) uses a different construction in two dimensions with explaining the 
behavior of the light rays but not the electromagnetic waves. The argument 
of [93] is only valid outside the cloaked region; it doesn't take into account 
the behavior of the waves on the entire region, including the cloaked re- 
gion and its boundary, the cloaking surface. In fact, the sequel [27], which 
gave numerical simulations of the electromagnetic waves in the presence of 
a cloak, states: "Whether perfect cloaking is achievable, even in theory, is 
also an open question" . In |35j we established that perfect cloaking is indeed 
mathematically possible at any fixed frequency. 

Before we discuss the paper [35] and other developments, we would like to 
point out that it is still an open question whether visual cloaking is feasi- 
ble in practice, i.e., whether one can realize such theoretical blueprints for 
cloaking over all, or some large portion of, the visible spectrum. The main 
experimental evidence has been at microwave frequencies [99] , with a lim- 
ited version at a visible frequency |104] . While significant progress has been 
made in the design and fabrication of metamaterials, including recently for 
visible light [78j 1102] , metamaterials are nevertheless very dispersive and one 
expects them to work only for a narrow range of frequencies. Even theoreti- 
cally, one can unfortunately not expect to actually cloak electromagnetically 
at all frequencies, since the group velocity cannot be faster than the velocity 
of light in a vacuum. 

In [351, Thm - 

was extended to the Helmholtz equation, which models 
scalar optics (and acoustic waves [231 [29] an d quantum waves under some con- 
ditions |125j ). and Maxwell's equations, corresponding to invisibility for gen- 
eral electromagnetic waves. The case of acoustic or electromagnetic sources 
inside and outside the cloaked region, leading to serious obstacles to cloaking 
for Maxwell's equations, was also treated. 



In Sec. 4.3, we consider acoustic cloaking, i.e., cloaking for the Helmholtz 
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equation at any non-zero frequency with an acoustic source p, 



(A g + k 2 )u = p, in 5(0, 2). (9) 

Physically, the anisotropic density is given by Ig^^g 1 -* and the bulk modulus 
byM 1 / 2 . 

For acoustic cloaking, even with acoustic sources inside 5(0, 1), we consider 
the same singular metric considered for electrostatics. However, we need to 
change the notion of a solution since for a generic frequency a H 1 (B (0,2)) 
smooth solution of the Helmholtz equation cannot simultaneously satisfy a 
homogeneous Neumann condition on the surface of the cloaked region 
Thm. 3.5] and have a Dirichlet boundary value that is a non-zero constant. 
We change the notion of solution for Helmholtz equation to a finite energy so- 
lution (see Sec. 4.3). The key ingredient of the rigorous justification of trans- 
formation optics is then a removable singularities theorem for the Laplacian 
on ^(5(0,2) \0). 



In Sec. 4.4 we consider the case of Maxwell's equations. In the absence of 
internal currents, the construction of [lUHS], called the single coating in [35J, 
still works once one makes an appropriate definition of finite energy solutions. 
However, cloaking using this construction fails in the presence of sources 
within the cloaked region, i.e., for cloaking of active objects, due to the 
nonexistence of finite energy, distributional solutions. This problem can be 
avoided by augmenting the external metamaterial layer with an appropriately 
matched internal one in D; this is called the double coating; see Sec. 4.4| 

In Sec. 5 we consider another type of transformation optics-based device, 
an electromagnetic wormhole. The idea is to create a secret connection be- 
tween two points in space so that only the incoming and the outgoing waves 
are visible. One tricks the electromagnetic waves to behave as though they 
were propagating on a handlebody, giving the impression that the topology 
of space has been changed. Moreover, one can manipulate the rays travel- 
ling inside the handle to obtain various additional optical effects; see Fig. 
3. Mathematically this is accomplished by using the single coating construc- 
tion with special boundary conditions on the cloaking surface. The main 
difference is that, instead of a point, we blow up a curve, which in dimension 
3 or higher is also an H l removable singularity for solutions of Maxwell's 
equations. 

Both the anisotropy and singularity of the cloaking devices present serious 
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Figure 3: An electromagnetic wormhole is obtained by blowing up a metric 
near a curve. This corresponds to e and /z on the exterior of a thickened 
cylinder causing electromagnetic waves to propagate as if a handle were at- 
tached to Euclidean space. Behavior of light rays: (Left) Rays travelling 
outside wormhole. (Right) A ray transiting wormhole. 
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challenges in trying to physically realize such theoretical plans using meta- 
materials. In Sec. 7, we give a general method, isotropic transformation 
optics, for dealing with both of these problems; we describe it in some detail 
in the context of cloaking, but it should be applicable to a wider range of 
transformation optics-based designs. 

A well known phenomenon in effective medium theory is that homogeniza- 
tion of isotropic material parameters may lead, in the small-scale limit, to 
anisotropic ones [ST]. Using ideas from [2J [2B] and elsewhere, we showed in 
[4*0! HU 22] how to exploit this to find cloaking material parameters that are 
at once both isotropic and nonsingular, at the price of replacing perfect cloak- 
ing with approximate cloaking of arbitrary accuracy. This method, starting 
with transformation optics-based designs and constructing approximations 
to them, first by nonsingular, but still anisotropic, material parameters, and 
then by nonsingular isotropic parameters, seems to be a very flexible tool 
for creating physically realistic designs, easier to implement than the ideal 
ones due to the relatively tame nature of the materials needed, yet essentially 
capturing the desired effect on waves for all practical purposes. 

In Sec. 8 we consider some further developments and open problems. 

2 Visibility for electrostatics: 
Calderon's problem 

Calderon's inverse conductivity problem, which forms the mathematical foun- 
dation of Electrical Impedance Tomography (EIT), is the question of whether 
an unknown conductivity distribution inside a domain in IR n , modelling, e.g., 
the Earth, a human thorax, or a manufactured part, can be determined from 
voltage and current measurements made on the boundary. A. P. Calderon's 
motivation for proposing this problem was geophysical prospection. In the 
1940's, before his distinguished career as a mathematician, Calderon was 
an engineer working for the Argentinian state oil company. Apparently, 
Calderon had already at that time formulated the problem that now bears 
his name, but he did not publicize this work until thirty years later [19]. 

One widely studied potential application of EIT is the early diagnosis of 
breast cancer [25] . The conductivity of a malignant breast tumor is typically 
0.2 mho, significantly higher than normal tissue, which has been typically 



10 



measured at 0.03 mho. See the book [4TJ] and the special issue of Physiological 
Measurement [51] for applications of EIT to medical imaging and other fields, 
and [13] for a review. 




Figure 4: Left: An EIT measurement configuration for imaging objects in 
a tank. The electrodes used for measurements are at the boundary of the 
tank, which is filled with a conductive liquid. Right: A reconstruction of 
the conductivity inside the tank obtained using boundary measurements. 
[Jari Kaipio, Univ. of Kuopio, Finland; by permission.] 

For isotropic conductivities this problem can be mathematically formulated 
as follows: Let Q be the measurement domain, and denote by o~(x) the 
coefficient, bounded from above and below by positive constants, describing 
the electrical conductivity in Q. In Q the voltage potential u satisfies a 
divergence form equation, 

V • aVu = 0. (10) 

To uniquely fix the solution u it is enough to give its value, /, on the 
boundary. In the idealized case, one measures, for all voltage distributions 
u\dn = f on the boundary the corresponding current fluxes, v- aVu, over the 
entire boundary, where v is the exterior unit normal to dQ. Mathematically 
this amounts to the knowledge of the Dirichlet-to-Neumann (DN) map, A a , 
corresponding to a, i.e., the map taking the Dirichlet boundary values of the 



solution to (10) to the corresponding Neumann boundary values, 
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A CT : u\ 9n i ^ v- aVu\m- (11) 
Calderon's inverse problem is then to reconstruct a from A a . 

In the following subsections, we give a brief overview of the positive results 
known for Calderon's problem and related inverse problems. 

A basic distinction, important for understanding cloaking, is between isotropic 
conductivities, which are scalar-valued, and anisotropic conductivities, which 
are symmetric matrix- or tensor-valued, modelling situations where the con- 
ductivity depends on both position and direction. Of course, an isotropic 
o~(x) can be considered as anisotropic by identifying it with a(x)I nxn . 

Unique determination of an isotropic conductivity from the DN map was 
shown in dimension n > 2 for C 2 conductivities in [108J. At the writing of 
the current paper this result has been extended to conductivities having | 
derivatives in [H] and |92j . In two dimensions the first unique identifiability 
result was proven in [85] for C 2 conductivities. This was improved to Lip- 
schitz conductivities in [15] and to merely L°° conductivities in [5]. All of 
these results use complex geometrical optics (CGO) solutions, the construc- 



tion of which we review in Sec. |2.1| We briefly discuss in Sec. |2.2| shielding, 
a less satisfactory variant of cloaking which is possible using highly singular 
isotropic materials. 



In Sec. 2J3 we discuss the case of anisotropic conductivities, i.e., conductiv- 
ities that may vary not only with location but also on the direction. In this 
case, the problem is invariant under changes of variables that are the identity 
at the boundary. We review the positive results that are known about the 
Calderon problem in this setting. The fact that the anisotropic conductivity 
equation is invariant under transformations plays a crucial role on the con- 
structions of electromagnetic parameters that make objects invisible, but for 
those one needs to make a final leap to using singular transformations. 



2.1 Complex geometrical optics solutions 



In this section, we consider isotropic conductivities. If u is a solution of ( 10 ) 
with boundary data /, the divergence theorem gives that 

QM)-= [ o-\Vu\ 2 dx= [ A a (f)fdS (12) 
Jn Jan 
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where dS denotes surface measure. In other words Q a (f) is the quadratic 
form associated to the linear map A a (f), i.e., to know A. a (f) or Q a {f) for 
all / G Hi(dQ) is equivalent. The form Q a {f) measures the energy needed 
to maintain the potential / at the boundary. Calderon's point of view in 
order to determine a in Q was to find enough solutions u G H l (Q) of the 
conductivity equation div(crVw) = so that the functions | Vw| 2 span a dense 
set (in an appropriate topology). Notice that the DN map (or Q a ) depends 
non-linearly on a. Calderon considered the linearized problem at a constant 
conductivity. A crucial ingredient in his approach is the use of the harmonic 
complex exponential solutions: 



u = e x ' p , where p G C n with p ■ p = 0. (13) 

Sylvester and Uhlmann |108j constructed in dimension n > 2 complex geo- 
metrical optics (CGO) solutions of the conductivity equation for C 2 conduc- 
tivities similar to Calderon's. This can be reduced to constructing solutions 
in the whole space (by extending a = 1 outside a large ball containing fl) 
for the Schrodinger equation with potential. We describe this more precisely 
below. 

Let a G C 2 (lR rt ), a strictly positive in W 1 and a = 1 for \x\ > R, for some 
R > 0. Let L„u = V ■ aVu. Then we have 



a~2L a \a-2vj = (A - q)v, (14) 

where 

«=-)U (15) 



Therefore, to construct solutions of L a u = in IR n it is enough to construct 



solutions of the Schrodinger equation (A — q)v = with q of the form (15). 
The next result proven in [108J states the existence of complex geometrical 
optics solutions for the Schrodinger equation associated to any bounded and 
compactly supported potential. 



Theorem 2.1 Let q G L°°(lR n ), n > 2, with q{x) = for \x\ > R > 0. Let 

— 1 < 5 < 0. There exists e(S) and such that for every p G C n satisfying 

p-p = 
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and 

llfl + bR 1/2 allr.oom^ + 1 < 



11(1 + 


X 


2\ 1/2 „ II , % , i 
) 9||i°°(R«) + 1 









i/iere exists a unique solution to 

(A - g)v = 

of the form 

v = e x -»(l + ip q (x,p)) (16) 

withip q (-,p) e Lj(R n ). Moreover ifj q (-,p) e H%(R n ) and, for < s < 1, there 
exists C = C(n, s,5) > such that 

N^P)lk<7-^7 (17) 



IPI 



Here 



L 2 s (R n ) = {f ] J(l + \x\ 2 ) s \f(x)\ 2 dx<oo} 



with the norm given by \\f\\ 2 L 2 — J(l + \x\ 2 ) 5 \f(x)\ 2 dx and H™(M. n ) denotes 

5 

the corresponding Sobolev space. Note that for large \p\ these solutions 
behave like Calderon's exponential solutions e x ' p . The equation for ip q is 
given by 

(A + 2p.V)W = g(l+V ? ). (18) 



The equation (18) is solved by constructing an inverse for (A + 2p ■ V) and 
solving the integral equation 

^ = (A + 2p-V)- 1 (g(l + ^)). (19) 

Lemma 2.2 Let -1 < 5 < 0, < s < 1. Let p e C n \ 0, p ■ p = 0. 

Lei / e L 2 +1 (IR ?1 ). T/ien t/iere exzsfe a unique solution u p 6 L 2 (IR n ) o/ t/ie 
eonatzon 

A p u p := (A + 2p ■ V)n p = /. (20) 
Moreover u p G HjiW 1 ) and 

||Mp||fl-|(R») < ^ZT S 

for < s < 2 and /or some constant C s $ > 0. 
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The integral equation (19) with Faddeev's Green kernel [M] can then be 



solved in L 2 (IR n ) for large \p\ since 

(I - (A + 2p ■ Vy'q)^ = (A + 2p • V)- 1 *? 

and || (A + 2p ■ V) _1 g|| L 2^ L 2 < -A for some C > where || • ||iJ-+l| denotes 
the operator norm between L 2 (IR n ) and L 2 (IR n ). We will not give details of 



the proof of Lemma 2.2 here. We refer to the papers |108[ 1107] 



If is not a Dirichlet eigenvalue for the Schrodinger equation we can also 
define the DN map 

du , 



on 



where u solves 



(A - q)u = 0; u\ dQ = f. 

Under some regularity assumptions, the DN map associated to the Schrodinger 
equation A — q determines in dimension n > 2 uniquely a bounded poten- 
tial; see |108j for the smooth case, [87] for L°°, and |20j for potentials in a 
Fefferman-Phong class. 

The two dimensional results of |85j,|15j, [5 J use similar CGO solutions and the 
d method in the complex frequency domain, introduced by Beals and Coif- 
man in [9] and generalized to higher dimensions in several articles |10j.[T].|86j. 

More general CGO solutions have been constructed in [52] of the form 

u = e r(w) (a + r), (21) 

where V0 • = 0, |V0| 2 = |V?/>| 2 and is a limiting Carleman weight 
(LCW). Moreover a is smooth and non-vanishing and ||r||^2(Q) = O(-), 
IMI-ffHrj) = ^(1)- Examples of LCW are the linear phase <p(x) = x ■ G 
S 171-1 , used in the results mentioned above, and the non-linear phase <f)(x) = 
ln\x — xq\, where Xq G R n \ ch (O) {ch(-) denoting the convex hull) which 
was used in [52] for the problem where the DN map is measured in parts of 
the boundary. For a characterization of all the LCW in M n , n > 2, see |31j . 
In two dimensions any harmonic function is a LCW |113j . 

Recently, Bukhgeim [16] used CGO solutions in two dimensions of the form 



(21 ) with cf) = z 2 or <p = z 2 (identifying IR 2 ~ C) to prove that any compactly 
supported potential q G L p ,p > 2, is uniquely determined by Cauchy data of 
the associated Schrodinger operator. 
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Other applications to inverse problems using the CGO solutions described 
above with a linear phase are: 

• Quantum scattering: It is shown in [SI] and [SSJ that in dimension 
n > 2 the scattering amplitude at a fixed energy determines uniquely a 
two body compactly supported potential. This result also follows from 
|108] (see for instance [111] . |112j ). Applications of CGO solutions to 
the 3-body problem were given in |114] . In two dimensions the result of 
[TB] implies unique determination of the potential from the scattering 
amplitude at fixed energy 

• Scalar optics: The DN map associated to the Helmholtz equation 
A + k 2 n 2 (x) with an isotropic index of refraction n determines uniquely 
a bounded index of refraction in dimension 3 or larger, see e.g. [8111108] . 

• Optical tomography in the diffusion approximation: In this case 
we have V • a(x)Vu — a a (x)u — iuu = in Q where u represents the 
density of photons, a(x) the diffusion coefficient, and a a the optical 
absorption. Using the result of [108] one can show in dimension three 
or higher that if to ^ one can recover both a and a a from the cor- 
responding DN map. If u — then one can recover one of the two 
parameters. 

• Electromagnetics: The DN map for isotropic Maxwell's equations 
determines uniquely the isotropic electric permittivity, magnetic per- 
meability and conductivity [91J. This system can in fact be reduced to 
an 8 x 8 Schrodinger system, A • / 8x8 — Q [91]. 

For further discussion and other applications of CGO with linear phase solu- 
tions, including inverse problems for the magnetic Schrodinger operator, see 

mu. 

2.2 Quantum Shielding 

In [33], also using CGO solutions, we proved uniqueness for the Calderon 
problem for Schrodinger operators having a more singular class of potentials, 
namely potentials conormal to submanifolds of M. n ,n > 3. These may be 
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more singular than the potentials in [20] and, for the case of a hypersurface 
S, can have any strength less than the delta function 5s- 

However, for much more singular potentials, there are counterexamples to 
uniqueness. We constructed a class of potentials that shield any information 
about the values of a potential on a region D contained in a domain Q from 
measurements of solutions at dQ. In other words, the boundary information 
obtained outside the shielded region is independent of q\jj. On Q\D, these 
potentials behave like q(x) ~ —Cd(x,dD)~ 2 ~ e where d denotes the distance 
to 3D and C is a positive constant. In D, Schrodinger's cat could live for- 
ever. From the point of view of quantum mechanics, q represents a potential 
barrier so steep that no tunneling can occur. From the point of view of optics 
and acoustics, no sound waves or electromagnetic waves will penetrate, or 
emanate from, D. However, this construction should be thought of as shield- 
ing, not cloaking, since the potential barrier that shields q\o from boundary 
observation is itself detectable. 



2.3 Anisotropic conductivities 

Anisotropic conductivities depend on direction. Muscle tissue in the human 
body is an important example of an anisotropic conductor. For instance 
cardiac muscle has a conductivity of 2.3 mho in the transverse direction and 
6.3 in the longitudinal direction. The conductivity in this case is represented 
by a positive definite, smooth, symmetric matrix a = (a %3 (x)) on Q. 

Under the assumption of no sources or sinks of current in Q, the potential u 
in Q, given a voltage potential / on dQ, solves the Dirichlet problem 



t,3=l v 7 



OonO 



(22) 



The DN map is defined by 



a du 



(23) 



where v = {y l , . . . , v n ) denotes the unit outer normal to dfl and u is the solu- 



tion of (22 ). The inverse problem is whether one can determine a by knowing 
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A CT . Unfortunately, A CT doesn't determine a uniquely. This observation is due 
to L. Tartar (see [M] for an account). 

Indeed, let ip : Q — > Q be a C°° diffeomorphism with ip\gn = Id, the identity 
map. We have 

Aa = K (24) 
where a = ip*a is the push-forward of conductivity a in ijj, 

Here Difj denotes the (matrix) differential of tp, (Dip) T its transpose and the 



composition in (25) is to be interpreted as multiplication of matrices. 

We have then a large number of conductivities with the same DN map: any 
change of variables of Q that leaves the boundary fixed gives rise to a new 
conductivity with the same electrostatic boundary measurements. 

The question is then whether this is the only obstruction to unique iden- 
tifiability of the conductivity. In two dimensions, this was proved for C 3 
conductivities by reducing the anisotropic problem to the isotropic one by 
using isothermal coordinates |106] and using Nachman's isotropic result [85J. 
The regularity was improved in [105J to Lipschitz conductivities using the 
techniques of P2] and to L°° conductivities in [B] using the results of [3]. 

In the case of dimension n > 3, as was pointed out in [72], this is a prob- 
lem of geometrical nature and makes sense for general compact Riemannian 
manifolds with boundary. 

Let (M , g) be a compact Riemannian manifold with boundary; the Laplace- 
Beltrami operator associated to the metric g is given in local coordinates by 
(IT]). Considering the Dirichlet problem ^ associated to ([T]), we defined in 
the introduction the DN map in this case by 



n Pi 

i,j=i J 

where v is the unit-outer normal. 

The inverse problem is to recover g from A g . 



(26) 

an 
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If if) is a C°° diffeomorphism of M which is the identity on the boundary, and 
ip* g denotes the pull back of the metric g by ip, we then have that Q holds. 

In the case that M is an open, bounded subset of M n with smooth boundary, 
it is easy to see ([72]) that for n > 3 

A g = A CT , (27) 

where 

gij = ^/{n^jj^ = \ g \\gij m (28) 

In the two dimensional case there is an additional obstruction since the 
Laplace-Beltrami operator is conformally invariant. More precisely we have 

1 

a9 ~ a 9 

for any function a, a > 0. Therefore we have that (for n = 2 only) 

K {r 9 ) = A 9 (29) 
for any smooth function a > so that a\gM = 1- 

Lassas and Uhlmann [69] proved that ^ is the only obstruction to unique 
identifiability of the conductivity for real-analytic manifolds in dimension n > 



3. In the two dimensional case they showed that (29) is the only obstruction 
to unique identifiability for C°°-smooth Riemannian surfaces. Moreover these 
results assume that A is measured only on an open subset of the boundary. 
We state the two basic results. 

Let T be an open subset of dM. We define for /, supp / CT 

A 9 ,r(/) = A g (/)|r. 



Theorem 2.3 (n > 3) Let {M,g) be a real-analytic compact, connected Rie- 
mannian manifold with boundary. Let T C dM be real-analytic and assume 
that g is real-analytic up to T. Then (A g ^,dM) determines uniquely (M,g). 



Theorem 2.4 (n = 2) Let (M,g) be a compact Riemannian surface with 
boundary. Let T C dM be an open subset. Then (A 9 r, dM) determines 
uniquely the conformal class of (M, g) . 
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Notice that these two results don't assume any condition on the topology of 
the manifold except for connectedness. An earlier result of [72] assumed that 



(M, g) was strongly convex and simply connected and T = dM. Theorem 2.3 



was extended in [70] to non-compact, connected real-analytic manifolds with 
boundary. The number of needed measurements for determination of the 
conformal class for generic Riemannian surfaces was reduced in [47J. It was 
recently shown that Einstein manifolds are uniquely determined up to isom- 
etry by the DN map |46j . 

In two dimensions the invariant form of the conductivity equation is given 
by 



dw g ((3V g )u := g- x l 2 d % (gWpjtdju) = (30) 

where (3 is the conductivity and div ff (resp. V s ) denotes divergence (resp. 
gradient) with respect to the Riemannian metric g. This includes the isotropic 
case considered by Calderon with g the Euclidian metric, and the anisotropic 
case by taking (g ij = -y ij and (5 = \g\ 1/2 ). It was shown in |105] for bounded 
domains of Euclidian space that the isometry class of (/?, g) is determined 
uniquely by the corresponding DN map. 

We remark that there is an extensive literature on a related inverse problem, 
the so-called Gelfand's problem, where one studies the inverse problem of 
determining a Riemannian manifold from the DN map associated to the 
Laplace-Beltrami operator for all frequencies, see [51] and the references cited 
there. 



3 Invisibility for Electrostatics 

The fact that the boundary measurements do not change, when a conduc- 
tivity is pushed forward by a smooth diffeomorphism leaving the boundary 
fixed, can already be considered as a weak form of invisibility. Different me- 
dia appear to be the same, and the apparent location of objects can change. 
However, this does not yet constitute real invisibility, as nothing has been 
hidden from view. 

In invisibility cloaking the aim is to hide an object inside a domain by sur- 
rounding it with (exotic) material so that even the presence of this object can 
not be detected by measurements on the domain's boundary. This means that 
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all boundary measurements for the domain with this cloaked object included 
would be the same as if the domain were filled with a homogeneous, isotropic 
material. Theoretical models for this have been found by applying diffeomor- 
phisms having singularities. These were first introduced in the framework of 
electrostatics, yielding counterexamples to the anisotropic Calderon problem 
in the form of singular, anisotropic conductivities in IR n , n > 3, indistin- 
guishable from a constant isotropic conductivity in that they have the same 
Dirichlet-to-Neumann map pBl H5] . The same construction was rediscovered 
for electromagnetism in [93], with the intention of actually building such a 
device with appropriately designed metamaterials; a modified version of this 
was then experimentally demonstrated in [99J. (See also [73J for a somewhat 
different approach to cloaking in the high frequency limit.) 

The first constructions in this direction were based on blowing up the me- 
tric around a point [70] . In this construction, let (M, g) be a compact 2- 
dimensional manifold with non-empty boundary, let Xo G M and consider 
the manifold 

M = M\{x } 

with the metric 

Um{X, Xq) 

where (Im(x,xo) is the distance between x and xq on (M,g). Then (M,g) 
is a complete, non-compact 2-dimensional Riemannian manifold with the 
boundary dM = dM. Essentially, the point Xq has been 'pulled to infinity" . 
On the manifolds M and M we consider the boundary value problems 

r , w f A su = in M, 

< „ Q , , and < u = f on dM, 

u— f on am, I J ~ ' 

L [ ueL°°{M). 

These boundary value problems are uniquely solvable and define the DN 
maps 

A-M,gf = d u u\ dM -, ^M t gf — dv u \dM 

where d u denotes the corresponding conormal derivatives. Since, in the two 
dimensional case functions which are harmonic with respect to the metric 
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g stay harmonic with respect to any metric which is conformal to g, one 
can see that Am, 9 = Aj^~. This can be seen using e.g. Brownian motion 
or capacity arguments. Thus, the boundary measurements for (M, g) and 
(M,7j) coincide. This gives a counter example for the inverse electrostatic 
problem on Riemannian surfaces - even the topology of possibly non-compact 
Riemannian surfaces can not be determined using boundary measurements 
(see Fig. 5). 

The above example can be thought as a "hole" in a Riemann surface that does 
not change the boundary measurements. Roughly speaking, mapping the 
manifold M smoothly to the set M \ B M (x , p), where B M (x , p) is a metric 
ball of M, and by putting an object in the obtained hole B M (x , p), one could 
hide it from detection at the boundary. This observation was used in (HI H5] , 
where "undetectability" results were introduced in three dimensions, using 
degenerations of Riemannian metrics, whose singular limits can be considered 
as coming directly from singular changes of variables. Thus, this construction 
can be considered as an extreme, or singular, version of the transformation 
optics of |116| . 

The degeneration of the metric (see Fig. 6), can be obtained by considering 
surfaces (or manifolds in the higher dimensional cases) with a thin "neck" 
that is pinched. At the limit the manifold contains a pocket about which 
the boundary measurements do not give any information. If the collapsing of 
the manifold is done in an appropriate way, we have, in the limit, a singular 
Riemannian manifold which is indistinguishable in boundary measurements 
from a flat surface. Then the conductivity which corresponds to this metric 
is also singular at the pinched points, cf. the first formula in (33). The 
electrostatic measurements on the boundary for this singular conductivity 
will be the same as for the original regular conductivity corresponding to the 
metric g. 

To give a precise, and concrete, realization of this idea, let B(0,R) C M. 3 
denote the open ball with center and radius R. We use in the sequel 
the set N = 5(0,2), the region at the boundary of which the electrostatic 
measurements will be made, decomposed into two parts, N± = 5(0,2) \ 
5(0, 1) and N 2 = 5(0, 1). We call the interface £ = dN 2 between JVi and 
N 2 the cloaking surface. 

We also use a "copy" of the ball 5(0,2), with the notation M 1 = 5(0,2), 
another ball M 2 = 5(0, 1), and the disjoint union M of M\ and M 2 . (We will 
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Figure 5: Blowing up a metric at a point, after [70]. The electrostatic boun- 
dary measurements on the boundary of the surfaces, one compact and the 
other noncompact but complete, coincide. 
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Figure 6: A typical member of a family of manifolds developing a singularity 
as the width of the neck connecting the two parts goes to zero. 



see the reason for distinguishing between N and M.) Let gjk = Sjk be the 
Euclidian metrics in Mi and M 2 and let 7 = 1 be the corresponding isotropic 
homogeneous conductivity. We define a singular transformation 

F 1 :M 1 \{0}^N 1 , F 1 (x) = (^ + 1)^- < \x\ < 2, (31) 

and a regular transformation (diffeomorphism) F 2 : M 2 1— > Af 2 , which for 
simplicity we take to be the identity map F 2 = Id. Considering the maps F\ 
and F 2 together, F = (F u F 2 ), we define a map F : M \ {0} = (M 1 \ {0}) U 
M 2 -> \ E. 

The push-forward g = -F*g of the metric g in M by F is the metric in N 
given by 

" Qpp Qpq 



dxi dx k 

p,q=l 



(32) 

x=F~ 1 (y) 

This metric gives rise to a conductivity a in N which is singular in Ni, 

~_ J \g\ 1/2 g jk iorxeNi, ( . 

° - \ V k for x G N 2 . [66) 
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Thus, F forms an invisibility construction that we call the "blowing up a 
point" . Denoting by (r, <ft, 9) i— > (r sin 9 cos <j), r sin 9 sin (f>, r cos #) the spheri- 
cal coordinates, we have 

/ 2(r - l) 2 sinfl \ 

= 2sinfl , 1 < \x\ < 2. (34) 

\ 2(sin#)- 1 / 

Note that the anisotropic conductivity a is singular degenerate on £ in the 
sense that it is not bounded from below by any positive multiple of I. (See 



[62] for a similar calculation.) The Euclidian conductivity <P in N 2 (33) could 
be replaced by any smooth conductivity bounded from below and above by 
positive constants. This would correspond to cloaking of a general object 
with non-homogeneous, anisotropic conductivity. Here, we use the Euclidian 
metric just for simplicity. 

Consider now the Cauchy data of all solutions in the Sobolev space H l (N) 
of the conductivity equation corresponding to a, that is, 

d(a) = {{u\ dN ,vaVu\ dN ) : ueH^N), V-aVw = 0}, 

where v is the Euclidian unit normal vector of dN. 

Theorem 3.1 flZJffl ) The Cauchy data of all H 1 -solutions for the conduc- 
tivities a and 7 on N coincide, that is, Ci(a) = 6*1(7). 

This means that all boundary measurements for the homogeneous conduc- 
tivity 7 = 1 and the degenerated conductivity a are the same. The result 
above was proven in |m US] for the case of dimension n > 3. The same basic 
construction works in the two dimensional case [62] • For a further study of 
the limits of visibility and invisibility in two dimensions, see [?]. 

Fig. 2 portrays an analytically obtained solution on a disc with conductivity 
a. As seen in the figure, no currents appear near the center of the disc, so 
that if the conductivity is changed near the center, the measurements on the 
boundary dN do not change. 

The above invisibility result is valid for a more general class of singular cloak- 
ing transformations, e.g., quadratic singular transformations for Maxwell's 
equations which were introduced first in [18J. A general class, sufficing at 
least for electrostatics, is given by the following result from [45J: 
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Theorem 3.2 Let Q C M. n , n > 3, and g = (g^) a smooth metric on Q 
bounded from above and below by positive constants. Let D CC Q be such 
there is a C°° -diffeomorphism F : Q\ {y} — > Q \ D satisfying F\qq = Id and 
such that 

dF(x) > CqI, det(dF(x)) > C\ dist Rn (35) 

where dF is the Jacobian matrix in Euclidian coordinates on W 1 and cq, c\ > 
0. Let g be a metric in Q which coincides with g = F*g in Q \ D and is an 
arbitrary regular positive definite metric in D mt . Finally, let o and a be the 



conductivities corresponding to g andTj, cf. (28). Then, 

Ci(a) = CiOt). 



The key to the proof of Thm. |3.2| is a removable singularities theorem that 
implies that solutions of the conductivity equation in Q \ D pull back by 
this singular transformation to solutions of the conductivity equation in the 
whole Q. 



Returning to the case f2 = iV and the conductivity given by (33), similar 
type of results are valid also for a more general class of solutions. Consider 
an unbounded quadratic form, A in L 2 (N, Ig^^dx), 



Az[u,v]= / aVu-Vvdx 
Jn 



defined for u,v e T>(A„) = C™(N). Let A„ be the closure of this quadratic 
form and say that 



V-crVw = in N 

is satisfied in the finite energy sense if there is Uq G ^(N) supported in Ni 
such that u — u G T>(A$) and 

Av[u — uq,v\ — — aVuo- Vvdx, for all v G V(A^). 
Jn 

Then the Cauchy data set of the finite energy solutions, denoted by 
Cf. e .{5) = |(w|ajv, V' &Vu\on) : u is a finite energy solution of V- a Vm = o| 
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coincides with the Cauchy data C/. e .(7) corresponding to the homogeneous 
conductivity 7 = 1, that is, 

C u Xa) = C f . e .(l). (36) 

This and analogous results for the corresponding equation in the non-zero 
frequency case, 

V • ffVii = Am, 

were considered in [35]. We will discuss them in more detail in the next 
section. 

We emphasize that the above results were obtained in dimensions n > 3. 
Kohn, Shen, Vogelius and Weinstein [62] have shown that the singular con- 
ductivity resulting from the same transformation also cloaks for electrostatics 
in two dimensions. 



4 Optical Invisibility: 

Cloaking at Positive Frequencies 

4.1 Developments in physics 

Two transformation optics-based invisibility cloaking constructions were pro- 
posed in 2006 [73l|93]. Both of these were expressed in the frequency domain, 
i.e., for monochromatic waves. Even though the mathematical models can 
be considered at any frequency, it is important to note that the custom de- 
signed metamaterials manufactured for physical implementation of these or 
similar designs are very dispersive; that is, the relevant material parameters 
(index of refraction, etc.) depend on the frequency. Thus, physical cloaking 
constructions with current technology are essentially monochromatic, work- 
ing over at best a very narrow range of frequencies. The many interesting 
issues in physics and engineering that this difficulty raises are beyond the 
scope of this article; see [75] for recent work in this area. 

Thus, we will also work in the frequency domain and will be interested in 
either scalar waves of the form U(x,t) = u(x)e lkt , with u satisfying the 
Helmholtz equation, 

(A + k 2 n 2 {x))u{x) = p(x), (37) 
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where p(x) represents any internal source present, or in time-harmonic elec- 
tric and magnetic fields E(x,t) = E(x)e tkt , H(x,t) = H(x)e lkt , with E, H 
satisfying Maxwell's equations, 

V x H = -ikeE + J, V x E = ikpH, (38) 

where J denotes any external current present. 

To review the ideas of [93] for electromagnetic cloaking construction, let 
us start with Maxwell's equations in three dimensions. We consider a ball 
5(0,2) with the homogeneous, isotropic material parameters, the permittiv- 
ity Eq = 1 and the permeability /i = 1. Note that, with respect to a smooth 
coordinate transformation, the permittivity and the permeability transform 



in the same way (25) as conductivity. Thus, pushing e and /xq forward 



by the "blowing up a point" map F\ introduced in (31) yields permittivity 
e{x) and permeability Jl(x) which are inhomogeneous and anisotropic. In 
spherical coordinates, the representations of e(x) and Jl(x) are identical to 



the conductivity a given in (34). They are are smooth and non-singular in 



the open domain N\ := 5(0,2) \ 5(0,1) but, as seen from (34), degener- 
ate as \x\ — > 1 + , i.e. at the cloaking surface S = {\x\ = 1}. One of the 
eigenvalues, namely the one associated with the radial direction, behaves as 
2(|x| 2 — l) 2 and tends to zero as \x\ — ► 1 + . This determines the electromag- 
netic parameters in the image of Fi, that is, in iVj.. In N 2 we can choose the 
electromagnetic parameters e(x),fi(x) to be any smooth, nonsingular ten- 
sors. The material parameters in N 2 correspond to an arbitrary object we 
want to hide from exterior measurements. 

In the following, we refer to N := N\ U iV2 U S = 5(0, 2) with the described 
material parameters as the cloaking device and denote the resulting specifica- 
tion of the material parameters on iV by e, /L As noted, the representations 



of e and JI on iV\ coincide with that of a given by (33) in spherical coor- 
dinates. Later, we will also describe the double coating construction, which 
corresponds to appropriately matched layers of metamaterials on both the 
outside and the inside of S. 



The construction above is what we call the single coating [35J. This the- 
oretical description of an invisibility device can, in principle, be physically 
realized by taking an arbitrary object in N 2 and surrounding it with spe- 
cial material, located in Ni, which implements the values of e, jl. Materials 
with customized values of e and fi (or other material parameters) are referred 
to as metamaterials, the study of which has undergone an explosive growth 
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in recent years. There is no universally accepted definition of metamateri- 
als, which seem to be in the "know it when you see it" category. However, 
the label usually attaches to macroscopic material structures having a man- 
made one-, two- or three-dimensional cellular architecture, and producing 
combinations of material parameters not available in nature (or even in con- 
ventional composite materials), due to resonances induced by the geometry 
of the cells |115[ 133] . Using metamaterial cells (or "atoms", as they are 
sometimes called), designed to resonate at the desired frequency, it is pos- 
sible to specify the permittivity and permeability tensors fairly arbitrarily 
at a given frequency, so that they may have very large, very small or even 
negative eigenvalues. The use of resonance phenomenon also explains why 
the material properties of metamaterials strongly depend on the frequency, 
and broadband metamaterials may not be possible. 



4.2 Physical justification of cloaking 

To understand the physical arguments describing the behavior of electromag- 
netic waves in the cloaking device, consider Maxwell's equations exclusively 
on the open annulus Ni and in the punctured ball Mi \ {0}. Between these 
domains, the transformation F\ : Mi \ {0} — > Ni is smooth. Assume that 
the electric field E and the magnetic field H in Mi \ {0} solve Maxwell's 
equations, 

V x H = -ike E, Vx£ = ikii H (39) 

with constant, isotropic Eq, liq. Considering E as a differential 1-form E(x) = 
Ei(x)dx 1 + E 2 (x)dx 2 + E 3 (x)dx 3 we define the push-forward of E by F\, 
denoted E = (F^E, in N u by 

3 3 3 

E(x) = J2Ej(x)dx j = (^(- DF_1 )i^) EkiF^i^dx^ x = F(x). 
j=i j=i k=i 

Similarly, for the magnetic field H we define H = (Fi)*H in iV\. Then E 
and H satisfy Maxwell's equations in N\, 

V x H = -ikeE, Vx£ = ikJIH, (40) 
where the material parameters in e, fx are defined in N\ by 

e = (Fi)*£ = ct, V- = (Fi)*/i = o. 
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Here a is given by (34). 



Thus, the solutions (E, H) in the open annulus N\ and solutions (E, H) in the 
punctured ball Mi \ {0} are in a one-to-one correspondence. If one compares 
just the solutions in these domains, without considering the behavior within 
the cloaked region N2 or any boundary condition on the cloaking surface E, 
the observations of the possible solutions of Maxwell's equations at dN = 
dB(0, 2) are unable to distinguish between the cloaking device N, with an 
object hidden from view in N 2 , and the empty space M. 

One can also consider the behavior of light rays, corresponding to the high 
frequency limits of solutions; see also [73], which proposed cloaking for scalar 
optics in M 2 . These are, mathematically speaking, the geodesies on the man- 
ifolds (Mi, g) and (Ni,g), see Fig. 7. One observes that almost all geodesies 
fi on N\ don't hit the cloaking surface S but go around the domain (N2,g) 
and have the same intrinsic lengths (i.e., travel times) as the corresponding 
geodesies Jl = F 1 ~ 1 ( y u) on (M^g). Thus, roughly speaking, almost all light 
rays sent into N\ from dN go around the "hole" N2, and reach dN in the 
same time as the corresponding rays on M. 

The cloaking effect was justified in [93] on the level of the chain rule for Fx, 
and in the sequels jOU [27] on the level of rays and numerical simulations, on 
Ni. We will see below that studying the behavior of the waves on the entire 
space, including in the cloaked region N 2 and at the cloaking surface S, is 
crucial to fully understanding cloaking and its limitations. 

A particular difficulty is that, due to the degeneracy of £ and Jl, the weighted 
L 2 space defined by the energy norm 

\\E\\ 2 1 +\\H\\ 2 1 = / (^E.Ek + Jl^HjH^dx (41) 

includes forms, which are not distributions, i.e., not in the dual of the vector 
fields having C£°(N) coefficients. Indeed, this class contains the forms with 
the radial component behaving like 0((r — l) _a ) in the domain r > 1, where 
1 < a < 3/2. The meaning of the Helmholtz or Maxwell's equations for such 
"waves" is problematic, and to treat cloaking rigorously, one should consider 
the boundary measurements (or scattering data) of finite energy waves which 
also satisfy Maxwell's equations in some reasonable weak sense. Analysis of 
cloaking from this more rigorous point of view was carried out in [35] , which 
forms the basis for much of the discussion here. 
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Figure 7: Left, light rays are shown in the Euclidian space 1R 3 and right, the 
same light rays are shown when a cloaking device (N, g) is located in the ball 
-8(0, 1). The metamaterial in which the light rays travel is not shown; the 
sphere is the cloaking surface dB(0, 1). On left, the light rays correspond to 
geodesies on (Mi \ {0},g) and on right, the geodesies on (Nx,]f)- The map 
Fx maps the geodesies on Mi (not passing through origin) to those of N\. 
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4.3 Cloaking for the Helmholtz equation 



Let us start with the cases of scalar optics or acoustics, governed in the case 



of isotropic media by the Helmholtz equation (37). In order to work with 
anisotropic media, we convert this to the Helmholtz equation with respect to 
a Riemannian metric g. Working in dimensions n > 3, we take advantage of 



the one-to-one correspondence (28 ) between (positive definite) conductivities 



and Riemannian metrics g. Let us consider the Helmholtz equation 

(A g + k 2 )u = p, (42) 

where A g is the Laplace-Beltrami operator associated with the Euclidian 
metric gij = dij. Under a smooth diffeomorphism F, the metric g pushes 
forward to a metric g = F*g, and then, for u = u o F, we have 

(A g + k 2 )u = p ^ (A g + k 2 )u = p, (43) 

where p = p o F. 

Next we consider the case when F(x) is not a smooth diffeomorphism, but 



the one introduced by (31 ), if x G Mi \ {0} and identity, if x e M 2 . 

Let / G L 2 (N,dx) be a function such that supp (/) n E = 0. We now give 
the precise definition of a finite energy solution for the Helmholtz equation. 

Definition 4.1 Let g be the Euclidian metric on M and g = F*g be the 

singular metric on N \ E. A measurable function u on N is a finite energy 
solution of the Dirichlet problem for the Helmholtz equation on N, 

(Ag + k 2 )u = f onN, (44) 
u\dN = h, 

if 

u G L 2 {N, \g\ 1/2 dx); (45) 
Z\N\neHl oc (N\Z,dx); (46) 

l \g\V 2 gV diUdjU dx < oo, (47) 
Jn\e 

u\dN = h; 



32 



and, for all ip G C°°(N) with ip\ dN = 0, 

-{Diujdji; + k 2 u^\g\ 1/2 ]dx = [ f{x)4){x)\g\ 1/2 dx (48) 

where D~u = \g\ 1 ' 2 'g i ^diU is defined as a Borel measure defining a distribution 
on N . 

Note that the inhomogeneity / is allowed to have two components, f\ and 
/ 2 , supported in the interiors of N ly N 2 , resp. The latter corresponds to an 
active object being rendered undetectable within the cloaked region. On the 
other hand, the former corresponds to an active object embedded within the 
metamaterial cloak itself, whose position apparently shifts in a predictable 
manner according to the transformation this phenomenon, which also 
holds for both spherical and cylindrical cloaking for Maxwell's equations, 
was later described and numerically modelled in the cylindrical setting, and 
termed the "mirage effect" [126J. 

Next we consider the relation between the finite energy solutions on N and 
the solutions on M. 

Theorem 4.2 (fifty) Letu = (ux,u 2 ) : {M 1 \{0})UM 2 ^Randu: iV\£ -> 
K. be measurable functions such that u = uo F . Let f = (f\, f 2 ) '■ (Mi\{0})U 
M 2 — > K and f : N \ S — ► R be L 2 functions, supported away from G Mi 
and T, C N such that f = jo F. At last, leth:dN^R,h: dM 1 -> R be 
such that h = h o F\. 

Then the following are equivalent: 

1. The function u, considered as a measurable function on N , is a finite 



energy solution to the Helmholtz equation (44 ) wi th inhomogeneity f 



and Dirichlet data h in the sense of Definition 4-1 
2. The function u satisfies 

{A g + k 2 ) Ul = h on M Xl Ul \ dMl = h, (49) 

and 

(A g + k 2 )u 2 = f 2 onM 2 , gi k v 3 d k u 2 \ dhh =b, (50) 

with 6 = 0. Here U\ denotes the continuous extension of u\ from Mi \ 
{0} to Mi. 
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Ito. 



Moreover, ifu solves (49) and (50) withb ^ 0, then the function u 
N \ £ — > K, considered as a measurable function on N , is not a finite energy 
solution to the Helmholtz equation. 



As mentioned in §|lj and detailed in (39], this result also describes a struc- 
ture which cloaks both passive objects and active sources for acoustic waves. 
Equivalent structures in the spherically symmetric case and with only cloak- 
ing of passive objects verified were considered later in 

The idea of the proof of Thm. 



4.2 is to consider F\ and Fi as coordinate 



transformations. As in formula (43), we see that if u is a finite energy solution 
of the Helmholtz equation (44) in N then u\ = u o _F\, defined in M\ \ {0}, 



satisfies the Helmholtz equation (49) on the set M\ \ {0}. Moreover, as the 
energy is invariant under a change of coordinates, one sees that m|mi\{o} is in 
the Sobolev space H l (Mi\{0}) . Since the point {0} has Hausdorff dimension 
less or equal the dimension of M. 3 minus two, the possible singularity of U\ at 
zero is removable (see e.g., [00]), that is, U\ has an extension to a function 



defined on the whole ball M\ so that the Helmholtz equation (49) is satisfied 
on the whole ball. 



Let us next discuss the appe aran ce of the Neumann boundary condition in 

l / 2 g tj diV is 



(50). Observe that in Def. 



4.1 



the Borel measure D~v 



\9\ 



absolutely continuous with respect to the Lebesgue measure for all functions 
v G Cq°(N). We can approximate the finite energy solution u of equation 



(44) with source /, supported away from S, by such functions. This yields 
that the measure of the cloaking surface satisfies D~m(E) = 0. Thus, using 

integration by parts, we see for arbitrary ip 6 C^(N) that 



= lim / [{Diu)djifj - k 2 uip\g\ 1/2 ]dx 

£ ^°+ J B(Q,l+e)\B(Q,l-e) 

= hm (7 -/ ) [u 3 {\g\ l ^d^dS{x), (51) 

£ ^ u + \JdB(0,l+e) JdB(0,l-e)J 

where dS is the Euclidian surface area. Changing coordinates by : 



dB(0, 1 + e) — > dB(0, 2e) in the first integral in (51 ) and letting e — > in the 
second integral, we see that 



= lim 



ij) dS 



du 
5 dv 



ij) dS, 



(52) 



s- 
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where ip = (Fi)*ip is a bounded function on Mi and u% is the solution of 
(49) in Mi, hence smooth near 0. Here §^L denotes the interior normal 
derivative. Thus, the first integral in (52) over the sphere of radius 2e goes to 
zero as e — > yielding that the last integral must vanish. As tjj is arbitrary, 
this implies that u satisfies the homogeneous boundary condition on the 
inside of the cloaking surface S. We point out that this Neumann boundary 
condition is a consequence of the fact that the coordinate transformation F 
is singular on the cloaking surface S. See also |61j for the planar case. 



4.4 Cloaking for Maxwell's equations 

In what follows, we treat Maxwell's equations in non-conducting and lossless 
media, that is, for which the conductivity vanishes and the components of e, /i 
are real valued. Although somewhat suspect (presently, metamaterials are 
quite lossy), these are standard assumptions in the physical literature. We 
point out that Ola, Paivarinta and Somersalo [91] have shown that cloak- 
ing is not possible for Maxwell's equations with non-degenerate isotropic, 
sufficiently smooth, electromagnetic parameters. 

We will use the invariant formulation of Maxwell's equations. To this end, 
consider a smooth compact oriented connected Riemannian 3-manifold M, 
DM 7^ 0, with a metric g, that we call the background metric. Clearly, in 
physical applications we take Met 3 with g being the Euclidean metric go. 
Time-harmonic Maxwell's equations on the manifold M are equations of the 
form 

curlE(x) =ikB(x), (53) 
curl H{x) = -ikD(x) + J. (54) 

Here the electric field E and the magnetic field H are 1-forms and the electric 
flux D and the magnetic flux B are 2-forms, and curl is the standard exterior 
differential d. The external current J is considered also as a 2-form. The 
above fields are related by the constitutive relations, 

D(x)=e(x)E(x), B(x) = n(x)H(x), (55) 

where e and fi are linear maps from 1-forms to 2-forms. Thus, in local 
coordinates on M, we denote 

E = Ej(x)dx j , D = D 1 (x)dx 2 A dx 3 + D 2 (x)dx 3 A dx 1 + D 3 (x)dx 1 A dx 2 , 
H = Hj(x)dx j , B = B l (x)dx 2 A dx 3 + B 2 (x)dx 3 A dx 1 + B 3 (x)dx 1 A dx 2 . 
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Using these notations, the constitutive relations take the form = e^ k E k 
and = Li jk H k . 

Note, that in the case of a homogeneous Euclidian space, where e = 1, Lio = 
1, the operators e and fi correspond to the standard Hodge star operator 
* : fi 1 (IR 3 ) — > fi 2 (IR 3 ) corresponding to the Euclidian metric (go)jk — $jk- On 
an arbitrary manifold (M, g) it is always possible to define the permittivity 
£ and permeability fx, to be the Hodge star operator corresponding to the 
metric g. Then, in local coordinates on M, 

This type of electromagnetic material parameters, which has the same trans- 
formation law, under the change of coordinates, as the conductivity, was 
studied in |6j 



To introduce the material parameters e(x) and fx{x) in the ball iV = B(0, 2) C 
M. 3 that make cloaking possible, we start with the singular map F\ given by 



(31). We then introduce the Euclidean metric on N 2 and the metric g = F*g 



in Ni. Finally, we define the singular permittivity and permeability in iV 



(33) 



using the transformation rules ( 56 ) which lead to the formulae analogous to 



~j k _ _ J Ig^^g^ for x e Ni, , . 

e ~ » ~ \ V k for x G N 2 . [b7) 

Clearly that, as in the case of Helmholtz equations, these material parameters 
are singular on E. 

We note that in N 2 one could define e and fx to be arbitrary smooth non- 
degenerate material parameters. For simplicity, we consider here only the 
homogeneous material in the cloaked region N 2 . 



4.5 Definition of solutions of Maxwell equations 

In the rest of this section, e = 1 and \x = 1 on the manifold M and e and fx 



are singular material parameter on iV defined in (57). 

Since the material parameters e and fx are again singular at the cloaking 
surface E, we need a careful formulation of the notion of a solution. 
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Definition 4.3 We say that (E,H) is a finite energy solution to Maxwell's 
equations on N , 

V x E = ik]l(x)H, V x H = -ike{x)E + J on N, (58) 

if E, H are one-forms and D := eE and B := J1H two-forms in N with 
L X {N ,dx)- coefficients satisfying 

[ e jk EjEkdVoix) < oo, (59) 
Jn 

I JP k H 3 ;H k dV {x) < oo; (60) 
Jn 

where dVo is the standard Euclidean volume and 

[ ((V x h) ■ E - ikh ■ $(x)H) dV (x) = 0, (61) 
Jn 

[ ((V x e) • H + e ■ (ike(x)E - J)) dV (x) = 
Jn 

for all 1- forms e,h on N having in the Euclidian coordinates components in 

c °°(iv). 



I FT II 2 



\H\ 



L2(jV,|g|VW (») 



Above, the inner product "•" denotes the Euclidean inner product. We em- 
phasize that in Def. |4.3 we assume that the components of the physical fields 
E,H,B, and D are integrable functions. This in particular implies that the 
components of these fields are distributions. Note that the map does not 
map distributions on M isomorphically to distributions on N. This is be- 
cause F* : (f> i-> (f> o F does not map Cq°(N) to Cq°(M). Hence, on M there 
are currents (i.e. sources) J, whose support contains the point zero that do 
not correspond to distributional sources J on N for which J = F* J in N\ E. 
Below we will show that in the case when a source J is not supported on 
N2 U S, there exist solutions for Maxwell's equations on iV with the corre- 
sponding source so that J = F* J in iV \ S. Also, we show that surprisingly, 
the finite energy solutions do not exist for generic currents J. Roughly speak- 
ing, the fact that the map F can not be extended to the whole M, so that it 
would map the differentiable structure on M to that of N, seems to be the 
reason for this phenomena. 

Below, we denote M \ {0} = (Mi \ {0}) U M 2 . 
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Theorem 4.4 (JEE/ ) Let E and H be 1 -forms with measurable coefficients on 
M \ {0} and E and H be 1-forms with measurable coefficients on N\H such 
that E = F*E, H = F*H . Let J and J be 2-forms with smooth coefficients 
on M \ {0} and N \ E ; that are supported away from {0} and S such that 
J = F,J. 

Then the following are equivalent: 

1. The 1-forms E and H on N satisfy Maxwell's equations 

Vx£ = ikJi(x)H, V x H = -ike{x)E + J on N, (62) 

V X E\qn = f 



in the sense of Definition 4-3\ 



2. The forms E and H satisfy Maxwell's equations on M, 

V x E = ikfi(x)H, V x H = -ike{x)E + J on M u (63) 
v x E\ dMl = f 

and 

V x E = ikft(x)H, V x H = -ike{x)E + J on M 2 (64) 
with Cauchy data 

v x E\ dM2 =b e , ux H\ 9M2 = b h (65) 
that satisfies b e = b = 0. 



Moreover, if E and H solve (65), (64), and (65) with non-zero b e orb , then 



the fields E and H are not solutions of Maxwell equations on N in the sense 
of Definition 

Let us briefly discuss the proof of this theorem. In Euclidian space, with 



e = 1 and fx = 1, Maxwell's equations (38) with J = and k ^ imply that 



the divergence of D and B fields are zero, or equivalently that 

V • (eE) = 0, V • (jiH) = 0. 
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Since e = 1 and // = 1, we obtain using (38) and the basic formulae of 
calculus, 



Jl2 

AE = Y^ tt^E = V(V ■£)-VxVx£ = 0-Vx (ik/iH) = -k 2 E. 



, dx] 
3 =i J 



This implies the Helmholtz equation (A + k 2 )E = 0. Thus, removable sin- 



gularity results similar to those used to prove Thm. 4.2 for the Helmholtz 



equation can be applied to Maxwell's equations to show that equations (62) 



on N imply Maxwell's equations (63) first on Mi \ {0} and then on all of 



Mi. Also, analogous computations to those presented after Thm. 4.2 for the 
finite energy solutions (E, H) of Maxwell's equations yield that the electric 
field E has to satisfy the boundary condition v x = on the inside 

of the cloaking surface. As E and H are in symmetric roles, it follows that 
also the magnetic field has to satisfy v x = 0. Summarizing, these 

considerations show that the finite energy solutions that are also solutions in 
the sense of distributions, have outside the cloaking surface a one-to-one cor- 
respondence to the solutions of Maxwell's equations with the homogeneous, 
isotropic e§ and \x§ on Mi, but inside the cloaking region must satisfy hidden 
boundary conditions at E~. 



Thm. |4.4| can be interpreted by saying that the cloaking of active objects is 
difficult since, with non-zero currents present within the region to be cloaked, 
the idealized model leads to non-existence of finite energy solutions. The 
theorem says that a finite energy solution must satisfy the hidden boundary 
conditions 

ux E = 0, ux H = on dN 2 . (66) 

Unfortunately, these conditions, which correspond physically to the so-called 
perfect electrical conductor (PEC) and perfect magnetic conductor (PMC) 
conditions, simultaneously, constitute an overdetermined set of boundary 
conditions for Maxwell's equations on N 2 (or, equivalently, on M 2 ). For 
cloaking passive objects, for which J = 0, they can be satisfied by fields 
which are identically zero in the cloaked region, but for generic J, including 
ones arbitrarily close to 0, there is no solution. The perfect, ideal cloaking 
devices in practice can only be approximated with a medium whose material 
parameters approximate the degenerate parameters e and /L For instance, 
one can consider metamaterials built up using periodic structures whose ef- 
fective material parameters approximate e'and /L Thus the question of when 
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the solutions exist in a reasonable sense is directly related to the question of 
which approximate cloaking devices can be built in practice. We note that if 



E and H solve (63), (64), and (65) with non-zero b e or b h , then the fields E 



and H can be considered as solutions to a set of non-homogeneous Maxwell 



equations on N in the sense of Definition 4.3 



V x E = ik[i(x)H + K sur f, V x H = —ike(x)E + J + J sur f on N, 

where K sur f and J sur f are magnetic and electric surface currents supported 
on E. The appearance of these currents has been discussed in [35l |371 1124] . 
We note that there are many possible choices for the currents J sur f and K sur j. 
If we include a PEC lining on E, that in physical terms means that we add 
a thin surface made of perfectly conducting material on E, the solution for 
the given boundary value / is the one for which the magnetic boundary 
current vanish, K sur f = and the electric boundary current J sur f is possibly 
non-zero. Introducing this lining on the cloaking surface E turns out to be 
a remedy for the non-existence results, and we will see that the invisibility 
cloaking then be allowed to function as desired. 

To define the boundary value problem corresponding to PEC lining, denote 
by CfP(iV) the space of functions / : N — > K such that f]^ and f\^ 2 are C°° 
smooth up to the boundary. 

Definition 4.5 We say that (E,H) is a finite energy solution to Maxwell's 
equations on N \ E with perfectly conducting cloaking surface, 

Vx£ = ikJi(x)H, V x H = -ike{x)E + J on N \ E, (67) 
v x E\v = 

if E, H are one-forms and D := eE and B := J1H two-forms in N with 



^(N^dx)- coefficients satisfying conditions (59 60), and equations (61) hold 
for all 1-forms e and h on N having in the Euclidian coordinates components 
in C%(N), vanishing near dN , and satisfing v x e|s = from both sides of 
E. 



With such lining of E, cloaking is possible with the following result, obtained 
similarly to Thm. 5 in [53] (cf. [351 Thm. 2 and 3]). 
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Theorem 4.6 Let E and H be 1 -forms with measurable coefficients on M\ 
{0} and E and H be 1- forms with measurable coefficients on N \ E such 
that E = F*E, H = F*H . Let J and J be 2-forms with smooth coefficients 
on M \ {0} and N \ E ; that are supported away from {0} and E such that 
J = F,J. 

Then the following are equivalent: 



1. The 1-forms E and H on N satisfy Maxwell's equations (61) in the 



sense of Definition 4-5 



2. The forms E and H satisfy Maxwell's equations on M , 

V x E = ik/i(x)H, V x H = -ike(x)E + J on M u (68) 

and 

V x E = ik/i(x)H, V x H = -ike{x)E + J on M 2 , (69) 
v x E\ 9 m 2 = 0. 



The above results show that if we are building an approximate cloaking device 
with metamaterials, effective constructions could be done in such a way that 
the material approximates a cloaking material with PEC (or PMC lining), 
which gives rise to the boundary condition on the inner part of E of the form 
v x E\g M2 — (or v x H\q M2 = 0). Another physically relevant lining is 
the so-called SHS (soft-and-hard surface) [53, EEJ SHI E] • Mathematically, 
it corresponds to a boundary condition on the inner part of E of the form 
E(X) = H(X) = 0, where X is a tangent vector field on E. It is particularly 
useful for the cloaking of a cylinder {(xi,x 2 ,x 3 ) G M 3 : (xi,x 2 ) G D}, D C 
IR 3 , when X is the vector J| in cylindrical coordinates, see [33], [31] • Further 
examples of mathematically possible boundary conditions on the inner part 
of E, for a different notion of solution, can be found in |119j . 

The importance of the SHS lining in the context of cylindrical cloaking is 
discussed in detail in [37]. In that case, adding a special physical surface 
on E improves significantly the behavior of approximate cloaking devices; 
without this kind of lining the fields blow up. Thus we suggest that the 
engineers building cloaking devices should consider first what kind of cloak 
with well-defined solutions they would like to approximate. Indeed, building 
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up a material where solutions behave nicely is probably easier than building 
a material with huge oscillations of the fields. 

As an alternative, one can avoid the above difficulties by modifying the basic 
construction by using a double coating. Mathematically, this corresponds to 
using an F = (F\, F 2 ) with both Fi, F 2 singular, which gives rise to a singular 
Riemannian metric which degenerates in the same way as one approaches E 
from both sides. Physically, the double coating construction corresponds 
to surrounding both the inner and outer surfaces of £ with appropriately 
matched metamaterials, see [35] for details. 

5 Electromagnetic wormholes 

We describe in this section another application of transformation optics which 
consists in "blowing" up a curve rather than a point. In [36, 38J a blueprint 
is given for a device that would function as an invisible tunnel, allowing 
electromagnetic waves to propagate from one region to another, with only 
the ends of the tunnel being visible. Such a device, making solutions of 
Maxwell's equations behave as if the topology of M 3 has been changed to 
that M 3 #(§ 2 xS 1 ), the connected sum of the Euclidian space R and the 
product manifold § 2 x S 1 . The connected sum is somewhat analogous to an 
Einstein-Rosen wormhole [32] in general relativity, and so we refer to this 
construction as an electromagnetic wormhole. 

We start by considering, as in Fig. 8, a 3-dimensional wormhole manifold, 
M = Mi U M 2 j ~, with components 

M 1 = R 3 \(B(0,1)UB(P,1)), 
M 2 = § 2 x [0,1]. 

Here ~ corresponds to a smooth identification, i.e., gluing, of the boundaries 
dMi and dM 2 . 

An optical device that acts as a wormhole for electromagnetic waves at a 
given frequency k can be constructed by starting with a two-dimensional 
finite cylinder 

T = S 1 x [0,L] C M 3 , 

taking its neighborhood K = {x G M 3 : dist(a;,T) < p}, where p > 
is small enough and defining N = M 3 \ K. Let us put the SHS lining on 
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Figure 8: A two dimensional schematic figure of wormhole construction by 
gluing surfaces. Note that the components of the artificial wormhole con- 
struction are three dimensional. 



the surface dK, corresponding to the angular vector field X = dg in the 
cylindrical coordinates (r, 9, z) in M 3 , and cover K with an invisibility cloak 
of the single coating type. This material has permittivity e and permeability 
ju described below, which are singular at OK. Finally, let 

U = {x: dist(x,K) > 1} C M 3 . 

The set U can be considered both as a subset of N, U C N((Z M 3 ) and of the 
introduced earlier abstract wormhole manifold M, U C M\. Let us consider 
the electromagnetic measurements done in U, that is, measuring fields E and 
H satisfying a radiation condition that corresponds to an arbitrary current J 
that is compactly supported in U . Then, as shown in |38j . all electromagnetic 
measurements in U C M and U G N coincide; that is, waves on the wormhole 
device (N, e, JZ) in M 3 behave as if they were propagating on the abstract 
wormhole manifold M. 

In Figures 3 and 9 we give ray-tracing simulations in and near the wormhole. 
The obstacle in Fig. 3 is K , and the metamaterial corresponding to e and Ji, 
through which the rays travel, is not shown. 

We now give a more precise description of an electromagnetic wormhole. Let 
us start by making two holes in M 3 , say by removing the open unit ball 
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Figure 9: Ray tracing simulations of views through the bores of two worm- 
holes. The distant ends are above an infinite chess board under a blue sky. 
On left, L « 1; on right, L ~ 1. Note that blue is used for clarity; the 
wormhole construction should be considered essentially monochromatic, for 
physical rather than mathematical reasons. 
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B\ = B(0, 1), and also the open ball B 2 = B(P, 1), where P = (0, 0, L) is a 
point on the z-axis with L > 3, so that BiC\B 2 = 0. The region so obtained, 
Mi = IR 3 \ (Bi U B2), equipped with the standard Euclidian metric go and 
a "cut" 7! = {(0,0, z) : 1 < z < L — 1}, is the first component Mi of the 
wormhole manifold. 

The second component of the wormhole manifold is a 3— dimensional cylin- 
der, M 2 = § 2 x [0, 1], with boundary dM 2 = (§ 2 x {0})U(§ 2 x {1}) := § 2 U§ 2 . 
We make a "cut" 72 = {^VP} x [0, 1], where NP denotes an arbitrary point 
in § 2 , say the North Pole. We initially equip M 2 with the product metric, 
but several variations on this basic design are possible, having somewhat 
different possible applications which will be mentioned below. 

Let us glue together the boundaries <9Mi and <9M 2 . The glueing is done so 
that we glue the point (0, 0, 1) G dB(0, 1) with the point NP x {0} and the 
point (0,0, L - 1) G dB(P, 1) with the point NP x {1}. Note that in this 
construction, 71 and 72 correspond to two nonhomotopic curves connecting 
(0, 0, l)~iVPx {0} to (0, 0, L - 1) ~ NP x {1}. Moreover, 7 = 7l U 72 will 
be a closed curve on M. 

Using cylindrical coordinates, (r, 9, z) 1— > (r cos 9, r sin 9, z), let N 2 = {(r, 6, z) : 
|r| < 1, z G [0,L]} n iV and N\ = N \ N 2 ; then consider singular transfor- 
mations Fj : Mj \ jj — > IR 3 , j = 1,2, whose images are Ni,N 2 , resp., see 
[38] for details. For instance, the map F\ can be chosen so that it keeps the 
^-coordinate the same and maps (z,r) coordinates by /1 : (z,r) — > (z',r'). 
In the Fig. 10 the map f\ is visualized. 

Together the maps F± and F 2 define a diffeomorphism F : M \ 7 —>■ N , 
that blows up near 7. We define the material parameters e and Jl on iV by 
setting e = F*e and /I = F*fi. These material parameters (having freedom 
in choosing the map F) give blueprints for how a wormhole device could be 
constructed in the physical space M. 3 . 

Possible applications of electromagnetic wormholes (with varying degrees of 
likelihood of realization!), when the metamaterials technology has sufficiently 
progressed, include invisible optical cables, 3D video displays, scopes for 
MRI-assisted medical procedures, and beam collimation. For the last two, 
one needs to modify the design by changing the metric g 2 on M2 = § 2 x [0, 1] . 
By flattening the metric on S 2 so that the antipodal point SP (the south 
pole) to NP has a neighborhood on which the metric is Euclidian, the axis 
of the tunnel N 2 will have a tubular neighborhood on which e, fi are constant 
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Figure 10: Above: A schematic figure of fx, representing F\, in the (z,r) 
plane. Its image P corresponds to Ni in (z, r) coordinates. Below: The 
sets Q and R correspond to N 2 and N. In the figure, R = Q U P which 
corresponds to iV = Ni U N 2 in IR 3 . 




Figure 11: Schematic figure. Left: Some rays enter the wormhole and come 
out from the other end so that they return near where the ray entered to 
the wormhole. Right: The corresponding ray in the complement iV of the 
obstacle K shown in the (z, r) coordinates. Note that there are also closed 
light rays. 
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isotropic and hence can be allowed to be empty space, allowing for passage 
of instruments. On the other hand, if we use a warped product metric on 
M 2 , corresponding to S> 2 x {z} having the metric of the sphere of radius r(z) 
for an appropriately chosen function r : [0, 1] — ► 1R + , then only rays that 
travel through N 2 almost parallel to the axis can pass all the way through, 
with others being returned to the end from which they entered. 



Remark 5.1 Along similar lines, we can produce another interesting class of 
devices, made possible with the use of metamaterials, which behave as if the 
topology of 1R 3 is altered. Let Mi = IR 3 \ -6(0, 1) endowed with the Euclidian 
metric g and M 2 be a copy of Mi. Let M be the manifold obtained by 
glueing the boundaries dM\ and dM 2 together. Then M can be considered 
as a C°° smooth manifold with Lipschitz smooth metric. Let N% = M 2 and 
F x : Mi -> Ni be the identity map, N 2 = 5(0,1) \ B(0,p) with p > 0, 
and finally N = R 3 \ B(0,p). Let F 2 : M 2 -> N 2 be the map F 2 (x) = 
(p+ (1 — p)|x| _1 )|x| _1 x. Together the maps F\ and F 2 define a map that can 
be extended to a Lipschitz smooth diffeomorphism F : M — > iV . As before, 
we define on N the metric g = F*g, and the permittivity e and permeability p 
according to formula (57). As N = IR 3 ^IR 3 , we can consider the (N, e, p) as a 
parallel universe device on which the electromagnetic waves on M 3 ^IR 3 can be 
simulated. It is particularly interesting to consider the high frequency case 
when the ray-tracing leads to physically interesting considerations. Light 
rays correspond to the locally shortest curves on N, so all rays emanating 
from Ni that do not hit dN\ tangentially then enter N 2 . Thus the light 
rays in N\ that hit dNi non-tangentially change the sheet N\ to N 2 . From 
the point of view of an observer in N\, the rays are absorbed by the device. 
Thus on the level of ray-tracing the device is a perfectly black body, or a 
perfect absorber. Similarly analyzing the quasi-classical solutions, the energy, 
corresponding to the non-tangential directions is absorbed, up to the first 
order of magnitude, by the device. Other, metamaterial-based constructions 
of a perfect absorber have been considered in (68]. We note that in our 
considerations the energy is not in reality absorbed as there is no dissipation 
in the device, and thus the energy is in fact trapped inside the device, which 
naturally causes difficulties in practical implementation. On the level of the 
ray tracing similar considerations using multiple sheets have been considered 
before in 1761. 



47 



6 A general framework: 

singular transformation optics 



Having seen how cloaking based on blowing up a point or blowing up a 
line can be rigorously analyzed, we now want to explore how more general 
optical devices can be described using the transformation rules satisfied by 
n, (p, \),e and p. This point of view has been advocated by J. Pendry and his 
collaborators, and given the name transformation optics |116j . As discussed 
earlier, under a nonsingular changes of variables F, there is a one-to-one 
correspondence between solutions u of the relevant equations for the trans- 
formed medium and solutions u = u o F of the original medium. However, 
when F is singular at some points, as is the case for cloaking and the worm- 
hole, we have shown how greater care needs to be taken, not just for the sake 
of mathematical rigour, but to improve the cloaking effect for more physi- 
cally realistic approximations to the ideal material parameters. Cloaking and 
the wormhole can be considered as merely starting points for what might be 
termed singular transformation optics, which, combined with the rapidly de- 
veloping technology of metamaterials, opens up entirely new possibilities for 
designing devices having novel effects on acoustic or electromagnetic wave 
propagation. Other singular transformation designs in 2D that rotate waves 
within the cloak [21], concentrate waves [96J or act as beam splitters [97] 
have been proposed. Analogies with phenomena in general relativity have 
been proposed in [71] as a source of inspiration for designs. 

We formulate a general approach to the precise description of the ideal mate- 
rial parameters in a singular transformation optics device, N C R 3 , and state 
a "metatheorem" , analogous to the results we have seen above, which should, 
in considerable generality, give an exact description of the electromagnetic 
waves propagating through such a device. However, we wish to stress that, as 
for cloaking [35] and the wormhole [3UI 138"] . actually proving this "result" in 
particular cases of interest, and determining the hidden boundary conditions, 
may be decidedly nontrivial. 

A general framework for considering ideal mathematical descriptions of such 
designs is as follows: Define a singular transformation optics (STO) design 
as a triplet (M.jMjJ 7 ), consisting of: 

(i) An STO manifold, M = (M, g, 7), where M = (Mi, . . . , M k ), the dis- 
joint union of n-dimensional Riemannian manifolds (Mj, gj), with or without 
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boundary, and (possibly empty) submanifolds jj C int Mj, with dim 7,, < 
ri-2; 

(ii) An STO device, M = (N,E), where N = \J k j= i N j C M n and E = 
U 3 -=i ^ji with Ej a (possibly empty) hypersurface in Nj] and 

(iii) A singular transformation T = (Fi, . . . , Fj,), with each 
Fj : Mj \ 7j — > iVj \ Ej a diffeomorphism. 

Note that N is then equipped with a singular Riemannian metric g, with 
<?| Nj = {Fj)*{9j), i n general degenerate on Ej. Reasonable conditions need 
to be placed on the Jacobians DFj as one approached jj so that the 7jj have 
the appropriate degeneracy, cf. [12J Thm.3]. 

In the context of the conductivity or Helmholtz equations, we can then com- 
pare solutions u on M. and u on A/", while for Maxwell we can compare fields 
(E, H) on Ai (with e and \i being the Hodge-star operators corresponding to 
the metric g) and (E, H) on M . For simplicity, below we refer to the fields 
as just u. 

Principle of Singular Transformation Optics, or "A Metatheorem 
about Metamaterials" : If (Ai,J\f, J 7 ) is an STO triplet, there is a 1-1 
correspondence, given by u = uo J 7 , i.e., u\uj — C^liv,-) ° Fj, between finite 
energy solutions u to the equation(s) on N ' , with source terms f supported 
on M \ E ; and finite energy solutions u on Ai, with source terms f = f ° T , 
satisfying certain "hidden" boundary conditions on dM = [L =1 dMj. 



7 Isotropic Transformation Optics 

The design of transformation optics (TO) devices, based on the transforma- 



tion rule (25), invariably leads to anisotropic material parameters. Further- 
more, in singular TO designs, such as cloaks, field rotators [21J, wormholes 
PSIEE], beam-splitters [97], or any of those arising from the considerations of 
the previous section, the material parameters are singular, with one or more 
eigenvalues going to or oo at some points. 

While raising interesting mathematical issues, such singular, anisotropic pa- 
rameters are difficult to physically implement. The area of metamaterials 
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is developing rapidly, but fabrication of highly anisotropic and (nearly) sin- 
gular materials at frequencies of interest will clearly remain a challenge for 
some time. Yet another constraint on the realization of theoretically per- 
fect ( or ideal in the physics nomenclature) TO designs is discretization: the 
metamaterial cells have positive diameter and any physical construction can 
represent at best a discrete sampling of the ideal parameters. 

There is a way around these difficulties. At the price of losing the theo- 
retically perfect effect on wave propagation that ideal TO designs provide, 
one can gain the decided advantages of being able to use discrete arrays of 
metamaterial cells with isotropic and nonsingular material parameters. The 
procedure used in going from the anisotropic, singular ideal parameters to 
the isotropic, nonsingular, discretized parameters involves techniques from 
the analysis of variational problems, homogenization and spectral theory. 
We refer to the resulting designs as arising from isotropic transformation op- 
tics. How this is carried out is sketched below in the context of cloaking; 
more details and applications can be found in jlOj |UJ H2] . 

The initial step is to truncate ideal cloaking material parameters, yielding a 
nonsingular, but still anisotropic, approximate cloak; similar constructions 
have been used previously in the analysis of cloaking [93 |371 |62j [21]. This 
approximate cloak is then itself approximated by nonsingular, isotropic pa- 
rameters. The first approximation is justified using the notions of T- and 
G-convergence from variational analysis (8[ [30] , while the second uses more 
recent ideas from [21 [3J, [26] . 

We start with the ideal spherical cloak for the acoustic wave equation. For 



technical reasons, we modify slightly the cloaking conductivity (33) by setting 
it equal to 25^ k on -8(0, 1), and relabel it as a for simplicity. Recall that a 
corresponds to a singular Riemannian metric g^ that is related to by 

a»(x) = |s(x)r/V(x), \g\ = (det[^']) 2 (70) 

where [g^ k {x)] is the inverse matrix of [gjk{x)} and \g{x)\ = det\gjk(x)]. The 
resulting Helmholtz equation, with a source term p, 

d ,, , , ,,/„ .., , d 



x) 



E if 

j,k=i 

u\dN = f, 



-yi — (\g( x )\Vi g i\ x )—u)+uj 2 u=p on iV, (71) 



can then be reinterpreted by thinking of a as a mass tensor (which indeed 
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has the same transformation law as conductivity under coordinate diffeo- 
morphisms ) and \g\z as a bulk modulus parameter; (71) then becomes an 
acoustic wave equation at frequency u with the new source pig] 1 / 2 , 



(v- crV + uj 2 \g\i j u = p(x)\g\* on N, 



(72) 



u\dN = /• 



This is the form of the acoustic wave equation considered in [231 EHl [39] . (See 
also [28J for d = 2, and |89j for cloaking with both mass and bulk modulus 
anisotropic.) To consider equation (72) rigorously, we assume that the source 



p is supported away from the surface S. Then the finite energy solutions u 



of equation (72) are defined analogously to Def. 4.1 Note that the function 
l l 2 appearing in (72) is bounded from above. 



\9 

Now truncate this ideal acoustic cloak: for each 1 < R < 2, let p 
and define F R :R 3 \ B(0, p) -> R 3 \ B(0, R) by 



2(R-1) 



x := F R (y) 



1 



\y\ 

2 



JL 



for \y\ > 2, 
for p < \y\ < 2. 



We define the corresponding approximate conductivity, <jr as 
jk/\ _ \ o- jk (x) for \x\ > R, 



a R (x) 



2S 3k , for Ixl < R, 



(73) 



where a^ k is the same as in the first formula in (33) or, in spherical coor- 
Note that then a^ k (x) 



dinates, (34) 



aorix) 



for \x\ > R, where 



o"o = 1 is the homogeneous, isotropic mass density tensor. Observe that, for 
each R > 1, a R is nonsingular, i.e., is bounded from above and below, but 
with the lower bound going to as R \ 1. Now define 



g R {x) = det {(?r{x)Y 



64|a;|- 4 (|x| - l) 4 for R < \x\ < 2, 
64, for |x| < R, 



cf. (|70|). Similar to (|72|), consider the solutions of 

l 

9i 
./'• 



(V- ctrV + u 2 g R /2 )u R 

u R\dN 



V2 • at 

g R ' p m N 



(74) 



(75) 
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As in Thm. 4.2 by considering Fr as a transformations of coordinates one 
sees that 



u R {x) 

,± 



v^(F^(x)), for R < \x\ < 2, 
Vr{x), for \y\ < R, 



with v R satisfying 



and 



(A + u 2 )v+{y) = p{F R {y)) in p < \y\ < 2, 

Vr\8B(0,2) = f, 



(v 2 + Au 2 )v- R { y ) = My), ™\y\<R- (76) 



Since <jr and gn are nonsingular everywhere, we have the standard transmis- 
sion conditions on Hr := {x : |x| = R}, 

ur\x r+ = u R \x R -, (77) 
e r - crR\7uR\x R+ = e r - o" fl Vu fl | Sfl _, (78) 

where e r is the radial unit vector and ± indicates when the trace on Hr is 
computed as the limit r — ► R ± . 

The resulting solutions, say for either no source, or for p supported at the 
origin, can be analyzed using spherical harmonics, and one can show that the 
waves v for the ideal cloak are the limits of the waves for the approximate 



cloaks, with the Neumann boundary condition in (50) for the ideal cloak 
emerging from the behavior of the waves v R for the truncated cloaks. This can 
be seen using spherical coordinates and observing that the trace of the radial 
component of conductivity from outside, <T] ? r '|s fl+ , goes to zero as R — > 1 
but the trace cx^ r |s H+ from inside stays bounded from below. Using this, we 



can see that the transmission condition ( 78 ) explains the appearance of the 
Neumann boundary condition on the inside of the cloaking surface. 

To consider general conductivities, we recall that for a conductivity 7- ?fc (x) 
that is bounded both from above and below, the solution of the boundary 



value problem (22) in N is the unique minimizer of the quadratic form 

Q^v) = / 7 Vw- Vvdx (79) 



N 
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over the functions v G H 1 ^) satisfying the boundary condition v\qn = f. 

We use the above to consider the truncated conductivities or. Note that at 
each point x E N the non-negative matrix <Jr(x) is a decreasing function 
of R. Thus the quadratic forms v i— > Q aR {v) are pointwise decreasing. As 
the minimizer v of the quadratic form Q aR {v) + (h,v)^ with the condition 
v |av = / is the solution of the equation 

V-a R Vv = h, v\ dN = f, 

we can use methods from variational analysis, in particular T-convergence 



e.g., pOj) to consider solutions of equation (75). Using that, it is possible 



to show that the solutions ur of the approximate equations ( 75 ) converge to 



the solution u of (72) for the general sources p not supported on E in the 
case when uo 2 is not an eigenvalue of the equation ( 72 ) . 



Next, we approximate the nonsingular but anisotropic conductivity or with 
isotropic tensors. One can show that there exist nonsingular, isotropic con- 
ductivities 7„ such that the solutions of 

M^r 1/2 V-7n(x)V + u 2 )u n =p onN, (80) 

U n \dN = /, 



tend to the solution of (75) as n — > oo. This is obtained by considering 
isotropic conductivities 7 n (x) = /i n (|x|) depending only on radial variable 
r = \x\, where h n oscillates between large and small values. Physically, this 
corresponds to layered spherical shells having high and low conductivities. 
As the oscillation of h n increases, these spherical shells approximate an ani- 
sotropic medium where the conductivity has much lower value in the radial 
direction than in the angular variables. Roughly speaking, currents can eas- 
ily flow in the angular directions on the highly conducting spherical shells, 
but the currents flowing in the radial direction must cross both the low and 
high conductivity shells. Rigorous analysis based on homogenization theory 
[21 EE] is used for ^ 2 e and one can see that, with appropriately chosen 
isotropic conductivities j n , the solutions u n converge to the limit ur. These 
considerations can be extended to all u 2 6 C \ D, where D C E_ is a dis- 
crete set, by spectral-theoretic methods, [55J. More details can be found in 

MM- 



Summarizing, considering equations (80) with appropriately chosen smooth 



isotropic conductivities 7 n and bulk moduli gu and letting n — > oo with 
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R = R(n) — > 1, we obtain Helmholtz equations with isotropic and non- 
singular mass and bulk modulus, whose solutions converge to the solution of 



the ideal invisibility cloak (72) 



A particularly interesting application of the above construction is to quantum 
mechanics. Zhang, et al. |125j described an anisotropic mass tensor m and 
a potential V which together act as a cloak for matter waves, i.e., solutions 
of the corresponding anisotropic Schrodinger equation. This ideal quantum 
cloak is the result of applying the same singular transform F as used for 
conductivity, Helmholtz and Maxwell, and applying it to the Schrodinger 
equation with mass tensor fho = 5^ k , Vq = 0. Due to the anisotropy of in, 
and the singularity of both in and V, physical realization would be quite chal- 
lenging. However, using the approximate acoustic cloak, one can describe an 
approximate quantum cloak that should be much easier to physically realize. 



An analogue of the reduction (14), (15) of the isotropic conductivity equa- 



tion to a Schrodinger equation can be carried out for the acoustic equation. 
Letting E = u 2 , ip n (x) = r ^j 2 (x)u n (x) , and 



V*(x) : = 7„- 1/2 V 2 7 y 2 (*) - E^gf + E, (81) 
one computes that ip n satisfies the Schrodinger equation 

(-A + Vf)ip n = Ei/> n in N. 
Furthermore, the family {V^} acts an approximate cloak at energy E: 



Theorem 7.1 Approximate quantum cloaking. Let W be a potential W G 
L°°(S(0, 1)), andE e R not be a Dirichlet eigenvalue of -A on N = B(0, 2), 
nor a Neumann eigenvalue of —A + W on 5(0,1). The DN operators at 
dN for the Schrodinger operators corresponding to the potentials W + 
converge to the DN operator corresponding to free space, that is, 

lim A w+V B(E)f = Ao{E)f 
in L 2 (dN) for any smooth f on dN. 

The convergence of the DN operators also implies convergence of the scatter- 
ing amplitudes fTity : \im n _ >O0 a w+V E(E,8',9) = cto(E,9',9). 
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Figure 12: Left: E not a Neumann eigenvalue; approximate quantum cloak. 
Matter wave passes almost unaltered. Right: E a Neumann eigenvalue; 
potential supports almost trapped state. 
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(Note that this is not a consequence of standard results from perturbation 
theory, since the do not tend to as n — ► oo. Rather, as n —> oo, the 
become highly oscillatory near £ and sup^, |V^f(a;)| — > oo as n — > oo.) 

On the other hand, when E is a Neumann eigenvalue of —A + on 5(0, 1), 
then V® supports almost trapped states, which correspond to matter waves 
(i.e., quantum mechanical particles) which reside in B(0, 1) with high prob- 
ability. See Fig. 12 and [41 J for more details and applications. 



Remark 7.2 Parameter distribution similar to (73) have been studied in the 
physics literature in the context of realistically achievable layouts of metama- 
terials approximating an ideal cloak. Using other, apparently only slightly 
different designs, one obtains in the limit other cloaking devices with enforced 
boundary conditions on the inside of E; see [ID], where approximate cloaks 
are specified which give rise instead to the Robin boundary condition. 



8 Further developments and open problems 

The literature on metamaterials, cloaking, and transformation optics is grow- 
ing rapidly . We briefly describe here only a few recent developments and 
remaining challenges. See [75] for a variety of perspectives. 

(a) Although the first description [4~4~| |4"5] of the cloaking phenomenon was 
in the context of electrostatics, no proposals of electrostatic metamaterials 
that might be used to physically implement these examples have been made 
to date. A proposal for metamaterials suitable for magnetostatics (cloak- 
ing for which is of course mathematically identical to electrostatics), and 
magnetism at very low frequencies, is in |120j . Since [99], there has been 
a push to obtain cloaking at higher frequencies, with the visual part of the 
electromagnetic spectrum an obvious goal. Progress has been reported in 
[TT1 11041 1781 1102] . However, broadband visual cloaking seems at this point to 
be far off. It should also be pointed out that serious skepticism concerning 
the practical advantages of transformation optics based cloaking over earlier 
techniques for reducing scattering has been expressed [5"9"] . 

(b) Other boundary conditions at the cloaking surface, analyzed in the 
time domain, based on Von Neumann's theory of self-adjoint extensions and 
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using a different notion of solution than that considered here, have been 
studied in [EE71 HE] EH] . See also [T2T] . 

(c) For simplicity, in cloaking we have mainly considered singular trans- 



formations which are affine linear in r. (See, however, Thm. 3.2 ) In situ- 
ations where the measurements are made further from cloaked object, [18] 
introduced, for spherical cloaking, transformations nonlinear in the radial 
variable in order to give better impedance matching with the surrounding 
media, and this was further explored for cylindrical cloaking in [122J. 

(d) Effective medium theory for metamaterials is in its early development, 
and seems to be particularly difficult for materials assembled from periodic 
or almost-periodic arrays of small cells whose properties are based on reso- 
nance effects. A physical (although mathematically nonrigorous) analysis of 
this kind of media is in |103j . which makes implicit assumptions about the 
smoothness of the fields which are violated when the fields experience the 
blow up demonstrated in [HE1 E7J • Some recent work on homogenization in 
this context is in E31. However, further efforts in this directions are needed. 



(e) Existing theories of cloaking deal predominantly with non-relativistic 
media; see, however, [73]. It seems that developing a theory compatible with 
the relativistic framework would be important. Similarly, transformation 
optics in the context of nonlinear media seems likely to become significant 
as metamaterial technology develops. 

(f) At the cloaking surface £ the cloaking metric g on N\ has a conical 
singularity in the sense of geometric scattering theory. It would be interesting 
to understand the relationship between cloaking and other transformation 
optics constructions on the one hand and geometric scattering on the other. 
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